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A lattice-gas model is suggested for describing the ordering phenomena in alkali-metal fullerides of face- 
centered-cubic structure assuming that the electric charge of alkali ions residing in either octahedral or 
tetrahedral sites is completely screened by the first-neighbor Ceo molecules. This approximation allows us to 
derive an effective ion-ion interaction. The van der Waals interaction between the ion and Ceo molecule is 
characterized by introducing an additional site energy at the tetrahedral sites. This model is investigated by 
using a three-sublattice mean-field approximation and a simple cluster- variation method. The analysis shows 
a large variety of phase diagrams when changing the site energy parameter. 
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I. INTRODUCTION 

Since the discovery of superconductivity in several 
alkali-metal fullerides a large variety of A2:C6o inter- 
calation compounds has been prepared and investigated 
(for recent reviews see papers [^,|| ) . In the face-centered- 
cubic (FCC) structure of pristine Cgo the molecules are 
bonded via van der Waals interaction. In this structure 
the tetrahedral and octahedral interstitial spaces are suf- 
ficiently large to accommodate alkali atoms without sig- 
nificant distortion of the lattice. The intercalated alkali 
atoms transfer their s electron to Cgo molecules, con- 
sequently, the lattice is contracted by the electrostatic 
interaction for small atoms (e.g. Na). 

Because of the charge transfer the electrostatic energy 
plays an important role in the formation of different or- 
dered structures, as well as in the distortion of the host 
lattice for large alkali concentration {x > 3). Fleming 
et al. Q have determined the electrostatic energies for 
RbjjCgo compounds assuming point charges for the Rb+ 
ions, and the contribution of the Cqq ion is expressed as 
the minimum Coulomb repulsion of x point charges on a 
sphere of radius R, with R taken to be the Cgo nucleon 
radius of 3.5 A. In the calculation of electrostatic ener- 
gies Rabe et al. Q have assumed a uniformly charged 
spherical shell with the same radius for the Gqq ion. 

In the present work we introduce a simple lattice-gas 
model to study the ordering phenomena in alkali ful- 
lerides. In this model the interstitial sites may be empty 
or occupied by one type of alkali ion. In agreement with 
the large size of fuUerenes, here we assume that an A+ al- 
kali ion is completely screened by uniformly distributing 
its transferred charge on the first neighbor Ceo molecules. 
For example, the alkali ion residing at a tetrahedral (oc- 
tahedral) site is surrounded by four (six) Cgo molecules 
with charges of — e/4 (— e/6). In other words, each inter- 



calated alkali atom transfers charges to its neighboring 
Cgo molecules independently of the position of the re- 
maining A"*". This simplification makes the derivation of 
an effective pair interaction possible, which is evidently 
a short range one. The difference between the sizes of 
tetrahedral and octahedral voids is considered by intro- 
ducing extra site energy related to the van der Waals 
interactions. This site energy parameter depends on the 
radius of the intercalated particle, consequently it is char- 
acteristic of the type of alkali atom. In the knowledge of 
the effective interactions (and site energy) characterized 
by a few parameters we are able to study the ordering 
within the framework of the lattice-gas (Ising) formal- 
ism. Accepting the spherical shell charge distribution on 
Ceo as suggested by Rabe et al. Q this model reproduces 
the ordered ground state energies given by them. A sim- 
ple mean-field analysis is carried out to demonstrate the 
richness of phase diagrams when varying two parameters 
(site energy and radius of the spherical shell). 

The above formalism can describe the ordered struc- 
tures observed in alkali fullerides with FCC structure. 
Namely, in pristine Ceo all the interstitial sites are empty, 
whereas in the superconducting compounds with nominal 
composition AaCgo both the octahedral and tetrahedral 
sites are occupied by alkali ions. Only the octahe- 

dral sites are occupied in ACeo forming NaCl structure 
for A=Rb and Cs. |^ Filling of the tetrahedral sites re- 
sults in the CaF2 structure observed by Rossiensky et al. 
|l0| in the Na^jCeo system. According to our analysis the 
above phases are found to be stable if one takes the to- 
tal energy of charged Ceo molecule into consideration as 
discussed later on. 

The early experimental investigations of the alkali in- 
tercalated fullerides are summarized by Zhu et al. [ pd| in a 
provisional phase diagram. Recently, Poirier and Weaver 
p2| and Kuzmany et al. p3| proposed a phase diagram 
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for the Ka;C6o system showing a eutectoid transformation 
from the homogeneous KCgo state to the composition of 
Ceo and K3C60 phases. 12 1 Our calculation confirms the 



existence of such a phase digram in a narrow range of 
parameters. 

The above model explains adequately the formation 
of different intercalated structures. Our purpose is not 
to argue that the present description gives the correct 
lattice-gas model for a given A^Cqq system but rather to 
emphasize some general features highlighted by this ap- 
proach because the model is adaptable for other metal- 
Ceo systems. By this means the present model provides 
a framework to classify these systems. The relevant fea- 
tures will be illuminated by displaying a set of distinct 
phase diagrams. 

In this paper we restrict ourselves to the systems cor- 
responding to Aa;C6o with rigid FCC structure. The 
present description consider the Ceo molecule as a spher- 
ical object, that is, the phenomena related to the orien- 
tational ordering are neglected. Both the interaction be- 
tween two Cgo molecules and the electronic structure are 
simplified therefore the present model can not describe 
the chain formation discovered very recently by Chauvet 
et al. g 

The lattice-gas model is derived in the subsequent sec- 



tion. In Sec. |I| we discuss the ground states suggested 
by the model for different values of site energy. A three- 
sublattice mean-field analysis is performed in Sec. IV to 
demonstrate the variety of phase diagrams. These re- 
sults are confirmed by a simple (triangular) version of 
the cluster-variation method (CVM) in Sec. Finally 
we summarize our conclusions and discuss the possible 
continuation of the present approach in Sec. 



II. THE MODEL 

The electronic structure calculations for alkali ful- 
lerides have indicated that the higher lying s bands of 
the metal are unoccupied and the electrons fill the empty 
tiu states of the Ceo molecule, so the alkali atoms are 
ionized and their valence electrons are accepted by the 
fuUerene molecule forming a negative ion. [iqjl^ In such 
a situation the Coulomb interaction has the most impor- 
tant contribution to the energy of the system. Several 
authors have already studied the effect of the Madelung 
energy on the charge-transfer solids binary [ p^jTol 
and pseudobinary alloys pO[ | and intercalated graphite 
compounds ]2l[ |. 

Statistical physics is able to give an adequate descrip- 
tion of the Ising systems with short range interaction but 
a great deal arises in handling the infinite range Coulomb 
interaction. An effective way of eliminating this problem 
is to construct a renormalized short range interaction. 
For this purpose one can, for example, use the cluster 
expansion of the Madelung energy. |19| In our treatment 
a well converging series of pair interactions is derived in 



order to avoid the difficulties arising from the long range 
order interactions among the ions. 

Calculations taking into account the polarizability of 
the fuUerene molecules in the host lattice ||2^ suggested 
that the energy of two charged Ceo molecules in the in- 
tercalated solid can be well described as the interaction 
between two point charges in a dielectric media charac- 
terized by e. 

According to that when calculating the Coulomb en- 
ergy of the system, the metal atoms and Ceo molecules 
are treated as point charges. The charge of the alkali 
atoms has been taken to be e and the charge of the Ceo 
has been strongly related to the occupation of the near- 
est neighbor interstitial sites. The charge distribution 
may easily be described in the lattice-gas formalism in- 
troducing the site occupation variables rji which is 1 if 
the i-th interstitial site is occupied and if this site is 
empty. Within this formalism the charge at site i is ev- 
idently qi — erji and the charge transferred to the a-th 
Ceo molecule can be written as: 



(1) 



where Na means the set of first neighbor tetrahedral and 
octahedral interstitial sites around the a-th Ceo molecule 
and Cj denotes the portion of charge transferred from 
the alkali atom residing at the j-th site to each near- 
est neighbor Ceo molecule. More precisely, Cj — 1/6 for 
octahedral and 1/4 for tetrahedral sites. Such a choice 
provides that each alkali ion is completely neutralized by 
its first neighbor Ceo molecules. The model assumes that 
the transferred charges accompanies the alkali ions when 
hopping from site to site. This charge distribution evi- 
dently satisfies the requirement of the neutrality of the 
crystal. 

The energy of the system can be given in the following 
form: 
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where EintiQa) is the intramolecular energy, the Greek 
and Latin indices run over the Ceo sites and interstitial 
sites of the FCC lattice respectively. Substituting the 
charges from (|l|) the Coulomb energy can be expressed 
as: 



Ec = 



a^f3 jeNa, keN^ "'^ 
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On exchanging the order of summations the energy can 
be separated into the contributions of onsite and pair 
interactions, namely 
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and A'^i denotes the set of first neighbor Cgo positions 
around interstitial site i. 

The magnitude of V^'j decreases rapidly with the dis- 
tance between sites i and j. One can realize in Eq. |^ the 
Coulomb interaction between two cluster formed by the 
i-th and j~th alkali ion as the central site and their near- 
est neighbor Cgo molecules. Since the clusters are neu- 
tral and the systems have cubic symmetry the monopole, 
dipole and the quadropole terms vanish. So there are 
only higher order multipole interactions between two 
sites. The good convergence properties of the pair inter- 
action is a consequence of our basic assumptions, viz. the 
nearest neighbor shell of the Cgo molecules completely 
screens out the alkali ion charge. 

In order to evaluate the intramolecular energy a series 
of UHF-MNDO calculations have been performed 
for up to a sixfold ionized single Cqq molecule where 
the characteristic bond lengths of the system were kept 
constant n = 1.400 Aand r2 = 1.440 A . 
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FIG. 1. Total energy of CgQ ion. Diamonds are the calcu- 
lated values, dashed line represents the fitted parabola. 

The total energy of the charged molecule as a function 
of the number of excess electrons can be fitted quite well 
with a parabola (see Fig. given by: 



Etot ^ A + Bx + Cx^ 



(7) 



where x is the number of extra electrons on the molecule 
and A = -7635. 02eV: B = -4.71926eV: C = 1.6262eV. 



Such an almost perfect parabolic behavior does not hold 
for other fuUerenes with lower symmetry. 

The quadratic energy dependence of the isolated 
molecule on the excess charges allows us to introduce 
an effective radius for the Cgo molecule having a value of 
K = 0.312a where the lattice constant a is given by exper- 
iment From previous calculations |2^j2j,^ R « 4.8 
A can be obtained. The differences between the present 
and the cited calculations are due to the different choice 
of the equilibrium bond lengths. Since the Ceo molecule 
is embedded in the crystal lattice the energy of the sin- 
gle molecule is modified by the surrounding media. Sup- 
posing dielectric screening mechanism one can apply the 
Onsager-Kirkwood theory of the electrostatic solvation 
p6| for estimating the effect of the lattice: 
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where tq is the radius of the cavity into which the 
molecule is embedded. Since the n dipole moment of 
the charged Ceo is zero choosing the radius of the cavity 
to be equal the effective radius the intramolecular energy 
can be written as: 



E,^t =A + Bci 



2eR 



(9) 



The first term gives a constant contribution to the to- 
tal energy and the linear term shifts the onsite energies. 
Both effects may be transformed out by rescaling energy. 
The spherical shell model proposed by Rabe et al. jj] 
with the original radius of the bucky-ball leads essen- 
tially different ordering properties. The quadratic term 
gives a T^j" contribution to the pair interaction between 
those alkali ions which transfer some charge to the same 
Ceo molecule, i.e.: 



(10) 



where a e A^iHTV. 
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Obviously, this contribution vanishes 
Furthermore, the 
increases the site energies too with 



when Tij exceeds a threshold value 



E 



(11) 



The contribution of the van der Waals interaction to 
the site energy is distinct for the tetrahedral and octahe- 
dral sites and depends on the size of the alkali ion. To 
characterize this contribution we introduce an additional 
site energy, i.e.: 



Eo for octahedral sites, 

et + SEt for tetrahedral sites 



(12) 



where SEt is considered as the difference between the 
two types of sites and is expected to be positive for large 
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alkali ions. In this expression the sum of and e" is 
indicated by to and e* for the octahedral and tetrahedral 
sites, respectively. The numerical calculations give to = 
-1.17817e2/£a and = -1.37843e2/ea (where e^/a « 1 
eV), i.e. the tetrahedral sites are primarily occupied for 
5Et = 0. 

We are now in a position to express the total energy of 
the system within the lattice-gas formalism. The follow- 
ing Hamiltonian defines the energy for any configuration 
of alkali ions described by the variables ry^: 

i? = ^ em, + V^jV^r|j , (13) 

i ij 

where Vij = ^i'j and e^, V/j, are defined above 

by Eqs. (H, (1) and ® . 

The effective pair interactions Vij vs. Tij/a are shown 
in Fig. ||. Vij has different values for the same distances; 
in other words, the pair interaction depends on the ori- 
entation of the pair for some values of Vij. Figure ^ 
demonstrates clearly that Vij becomes extremely weak 
if Vij > a . The attractive contributions of the effective 
pair interaction (for a/2 < Tij < a) indicates the stability 
of alkali intercalated fuUerides at large alkali concentra- 
tions. 

As a comparison V-j is also plotted in Fig. ||. The 
striking difference between Vij and V-j indicates the im- 
portant role of the electronic energy of Cg,^ ions as will 
be investigated in detail in the subsequent sections. 
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FIG. 2. Effective pair interaction as a function of ion-ion 
distance. Open diamonds represent the pure coulombic in- 
teraction (V^'j), closed diamonds refer to full interaction (Vij) 
including the energy of charged Ceo molecule. 

III. GROUND STATES 



the sublattices are labelled for octahedral and 1 or 2 for 
tetrahedral sites. In the three-sublattice formalism the 
states are characterized by a vector consisting of sub- 
lattice occupations. For example, the state (cto, cti, (72) 
denotes an alkali distribution where the sites of sublat- 
tice V are occupied with a probability a^. In general 
< cTi/ < 1, however, a^, =0 or 1 for ordered structures. 
These ordered states are well known in crystallography. 
For example, the states (1,0,0), (0,1,1) and (1,1,1) are 
equivalent to the NaCl, CaF2, and BiFa structures. 




FIG. 3. Arrangement of octahedral (medium size spheres) 
and the two tetrahedral (dark and light small spheres) inter- 
stitial sites in an FCC lattice formed by Ceo molecules (large 
spheres) . The triangle represents the basic cluster of our C VM 
approximation. 

One can easily see that according to the present model 
the Ceo molecules are uniformly charged for the ordered 
states defined above. Consequently, the electrostatic en- 
ergies are equivalent to the Madelung energies deter- 
mined previously. |^J^ In agreement with the expecta- 
tion the energies of ordered states are equivalent to those 
predicted by Rabe et al. when choosing R = 3.5 A. In 
the above model, however, we have chosen larger R as 
detailed in the previous section. This correction modifies 
the total energy of the ordered states as plotted in Fig. ^ 
for 5Et = 1.1 . As a result the present model suggests 
three stable ground states [(0,0,0), (1,0,0) and (1,1,1)], 
which are connected by solid lines in the figure. That 
is, the system with a nominal alkali content a; = 1 has 
(1,0,0) phase in contrast to the former model predicting 
a stable (0,1,0) phase. 



The electrostatic energies of some ordered structures 
have already been investigated. Now we concen- 

trate on analysing such FCC structures which may be de- 
scribed by introducing three sublattices as demonstrated 
in Fig. ||. The sites of each sublattice form a FCC struc- 
ture equivalent to the cage lattice. For later convenience 
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FIG. 4. Smallest electrostatic energies of the diflterent 
ordered structures of A^Ceo at fixed concentrations for 
SEt — l.le^/ea. The solid line joins the stable phases. 

In order to study the stability of the ground states 
we introduce the chemical potential /x which controls the 
alkali content (x) in the system. The stable state as a 
function of /i is determined by the minimum of Gibbs 
potential per Ceo defined as 

00 = Ho- MX (14) 

where the index refers to zero temperature and Ti.o 
denotes the total energy per Cgo- A simple numerical 
calculation gives the stable state as a function of fi and 
SEt- The results are summarized in Fig. ^ This map 
shows that the alkali content increases with the chemical 
potential for a fixed 6Et . The energy relations are similar 
to those plotted in Fig. | if SEt > l.OlOGe^/ea. There is 
a region of SEt when the (0,0,0) state transforms directly 
to (1,1,1) for increasing fj, as is suggested theoretically by 
Fleming et al. |§ If SEt < 0.1406eVea then the (0,1,1) 
state is stable. This behavior is a consequence of the 
fact that SEt does not give an energy contribution to the 
states (0,0,0) and (1,0,0) however, it shifts the energy of 
the states (0,1,1) and (1,1,1) by the same value. 

We have evaluated the energies of several ordered 
states related to a different sublattice division of the 
tetrahedral sites. In order to determine the ground 
states of the systems both the octahedral and tetrahedral 
(FCC) lattices have been decomposed into four simple 
cubic lattices and the energies of all the possible ordered 
structures on the 12 sublattices have been calculated. 

In comparison with previous results we obtained higher 
energies for all the states distinguishable from the former 
ones. This investigation confirms the above sublattice 
division thereby providing a framework for later analysis. 

It is mentioned here that the ground state energies are 
independent of the screening mechanism. The role of 
this mechanism becomes important when studying the 
sublattice occupations as a function of temperature. 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 
5E 

FIG. 5. Ground states of A^Ceo as a function of chemical 
potential fi and 5Et. 

IV. MEAN-FIELD APPROXIMATION 

In this section the thermodynamic properties of the 
lattice-gas model defined by the Hamiltonian ([l3| ) are 
investigated by using a three-sublattice mean-field ap- 
proximation. It is assumed that the interstitial points 
are occupied by an alkali ion with the same probability 
CT^ within the sublattice = 0, 1, 2. 

By means of the mean-field approximation the energy 
per Ceo molecule may easily be expressed by the variables 
ai, as 

7i = ^ e^a^, + ^ X! Ji^rCTyCFr (15) 

where the Jin coupling constants contain all the interac- 
tions between an alkali ion in sublattice v and those ions 
residing in sublattice r, i.e., 

where Sy denotes the set of sites belonging to sublattice 
V. The tensor of coupling constants is symmetric, that 
is, Jvr = Jrv Furthermore, Jqi = J02 and Jn — J22 
because of the equivalence of the tetrahedral sublattices. 
Consequently, we have only four different coupling con- 
stants characteristic of the model in mean-field approxi- 
mation. Their values may be determined by straightfor- 
ward numerical calculation which gives 

e2 

Joo = -1.42879— , 

ea 

Jii = -1.60386— , 

ea 

e2 

Joi = -0.29000— , 

ea 

J12 = -0.86560— (17) 

ea 
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if R/a = 0.312. Notice that all the J^^ are negative, i.e., 
the attractive forces rule over the short range repulsion 
(see Fig. I). . 

Here it is worth mentioning that Eq. permits some 
modification of pair interactions Vij leaving the mean- 
field parameters (as well as the ground state energy) 
unchanged. For example, one can remove the interac- 
tions for large distances adding suitable corrections to 
the short range terms. In some sense the above screening 
mechanism may be considered as an example of such a 
process. This simplification may be very useful for Monte 
Carlo simulations. As a limit case, one can introduce a 
model characterized by the first four Vij terms with a 
suitable choice of their strength. 

It is evident that the mean-field Hamiltonian ( |l5|) re- 
produces the energies of the ordered states discussed 
above. However, the coupling constants J^^ cannot be 
determined in the knowledge of the Madelung energies 
for all the ordered structure because this Hamiltonian 
contains more information, viz. it gives the potential en- 
ergies for the empty sites too. 

The thermodynamic properties of the system are de- 
termined by the minimization of Gibbs potential 



Q — H — TS — fi ^ a„ 



(18) 



with respect to ctj/ for fixed temperature T and chemical 
potential /i. In the above expression the configurational 
entropy per Cgo is given as 

S = -ks ^[(7, log a, + (1 - a,) log(l - a,)] (19) 



where fc^ is the Boltzmann constant. 

The results of the numerical calculations are summa- 
rized in temperature-composition (T — x) phase diagrams 
for different values of SEt . The alkali concentration is re- 
lated to sublattice occupations (x = '^cr,^) determined 
directly by the numerical procedure. 

In agreement with the ground state analysis. Fig. ^ il- 
lustrates that only the (0,0,0), (0,1,1) and (1,1,1) phases 
are stable at low temperatures for 6Et ~ 0. For interme- 
diate compositions the system segregates into two phases 
as indicated in the figure. It is emphasized that in the 
high temperature phase the tetrahedral and octahedral 
sites are occupied with different probabilities. More pre- 
cisely, the tetrahedral sites are preferred to the octahe- 
dral ones for SEt — 0. For example, if the temperature 
is decreased for a fixed concentration x = 2 the present 
model suggests a continuous ordering process character- 
istic of the CaF2-type superionic conductors. [M 
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FIG. 6. Phase diagram for SEt = in mean-field approxi- 
mation. 

The (0,1,1) structure becomes unstable at T = when 
SEt > 0.1406e^/£a. Fig. |^ shows a situation when the 
homogeneous (0,1,1) phase may be observed in the sys- 
tem but it decomposes into two states [(0,0,0)-f (1,1,1)] 
at low temperatures. A eutectoid transition appears at 
ksT/ie'^/ea) = 0.2442 for SEt = 0.15e^/ea. To demon- 
strate this eutectiod transition the value of SEt is chosen 
to be close to the threshold value of its stability at T = 0. 

The decomposition of the high temperature phase into 
the phases (0,0,0) and (1,1,1) may be observed in a region 
of SEt defined in Sec. 



Ill 



A typical phase diagram that 
characterizes the phase separation at finite temperatures 
will be shown in the subsequent section. 
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FIG. 7. Phase diagram in mean-field approximation for 
SEt = O.lSe^/ea exhibits a eutectoid phase transition. 

As mentioned above the energy of both the (0,1,1) 
and (1,1,1) phases increases with SEt meanwhile the en- 
ergy of the phases (0,0,0) and (1,0,0) remains unchanged, 
and the (1,0,0) phase appears in the phase diagrams 
for SEt > 1.0106e^/£a. More precisely, close to this 
stability threshold value one can observe a eutectoid 
transition (see Fig. ||) which is analogous to the situ- 
ation discussed above. This phase diagram is similar 
to those suggested by Poiriet and Weaver |l^ and Kuz- 
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many et al. for the K2;C6o system for a; < 3. With 
X-ray photoemission and Raman spectroscopy they ob- 
served a reversible transformation from KCgo to a mixed 
phase of Cgo and K3C60 at the eutectoid temperature 
{Tent = 150 ± 10 °C). For 5Et = e^/ea the present ap- 
proach predicts kBTeut — 0.250e^/£a. 

A new type of phase diagram will be characteristic of 
the system when the (1,0,0) phase becomes stable at zero 
temperature. Instead of displaying now a typical plot, 
m Sec. we compaire the phase diagrams suggested by 
mean-field approximation and CVM for 5Et = l.le^/ea. 



ous analysis of all the possible diagrams as a function of 
R and 5Et, however, goes beyond the scope of the present 
paper. 
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FIG. 8. Mean-field analysis of lattice-gas model suggests 
eutectoid phase transition for 5Et = e^/ea. 

It is emphasized that in the above phase diagrams (in- 
cluding those plotted in Sec. |^ both tetrahedral sub- 
lattices are occupied with the same probability implying 
the union of these sublattices. However, some modifica- 
tion of the model parameters results in the appearance 
of the (0,1,0) and (1,1,0) states. We have also evaluated 
some phase diagrams accepting the spherical shell model 
suggested by Rabe et al. with a radius of the Ceo 
molecule. In order to demonstrate the difference between 
these choices we display two phase diagrams correspond- 
ing to i? = 3.5 A. In the first example the (0,1,0) state 
is stable at T = as illustrated in Fig. ||. The addi- 
tional feature observable in this plot is the appearance of 
an order-disorder transition confined to tetrahedral sites 
for X = 1 meanwhile the octahedral sites are practically 
empty. This process is analogous to an "antiferromag- 
netic" to "paramagnetic" transition, in other words the 
(0,1,0) state transforms into (0,1/2,1/2) at a critical tem- 
perature (/cbT/v = 0.555e^/ea). 

The mean-field analysis suggests a very interesting 
phase diagram for 6Et = e^/a (see Fig. |l^). In this case, 
if the temperature for a; = 2 is decreased the high tem- 
perature phase decomposes into two phases whose alkali 
concentrations are close to 1 and 3. At lower tempera- 
tures, however, the (1,1,0) phase becomes stable in com- 
parison to the coexistence of states (1,0,0) and (1,1,1). 

The richness of the predicted phase diagrams is surpris- 
ing because we have varied only two parameters. Rigor- 




FIG. 9. Phase diagram for SEt — when the electronic 
energy of Cqq is approximated by a spherical shell model 
with radius 3.5 A. 
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V. CLUSTER- VARIATION METHOD 

The cluster variation methods (CVM) are successfully 
used for deriving phase diagrams with sufficient accu- 
racy. 1^ Two basic difficulties emerge when adapting 
this technique to the present lattice-gas model. On the 
one hand, the interstitial sites form a non-Bravais lat- 
tice; on the other hand, instead of the usual first (and 
second) neighbor interaction(s) there are a lot of param- 
eters Vij characteristic of the interaction between sites i 
and j. Here we suggest the application of a simple version 
of CVM to increase the accuracy of the three-sublattice 
mean-field approximation. 

First we introduce a set of variables i(??o, 'yi, '72) to 
denote the probability of a configuration (?7o, J?!, %) on 
three neighboring sites where r/^ refers to the occupa- 
tion of the site belonging to sublattice ly. The position 
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of such a three-point cluster is indicated by a triangle in 
Fig. ^. The probability of a configuration on a part of 
this triangle may be easily expressed by the values of t. 
For example, the probabilities for pair configurations are 
given as 



Poi{Vo,Vi) = ^t{vo,Vi,V2) 

m 

m 



(20) 



no 



and for single site configurations 

'11^,112 

S2{v2) = Y t{vo,m,m) 



(21) 



whose values are directly related to the sublattice occu- 
pations (Ty introduced previously, namely Si/(1) = cr^ and 

S^{0) = 1- Uy. 

The contribution of the first and second neighbor inter- 
actions to the system energy may be directly expressed 
by the variables Pi/r(l, 1) and only the remaining terms 
are approximated on the basis of mean-field theory. That 
is, now the energy per Cgo molecule is given as 

^(CVM) ^ ^v{\)[pm{l, 1) + P02(l, 1)] + 6F(2)pi2(l, 1) 

+ \Y.JlrS.{^)Sr{l)+Y.'uS.^^) (22) 

where ^(1) = 0.27459eVea and V{2) = -0.00386e7ea 
denotes the interactions Vij at the shortest and second 
shortest distances and J^^ is equivalent to J^n defined by 
Eq. (|l^) except the coupling constants should be reduced 
by the contributions of ViX) and ^(2), 



J'12 - 'h2 - 61/(2) 



(23) 



Here we restrict ourselves to demonstrating how the 
entropy may be derived on the analogy of pair approxi- 
mation developed by Bethe. ||2^ The non-Bravais lattice 
of interstitial sites may be built up by repeating period- 
ically the triangle of an octahedral and two tetrahedral 
sites as prescribed by the primitive cell of FCC structure. 
This construction makes it clear that each triangle is sur- 
rounded by twelve triangles with the same orientation. 
According to the pair approximation the entropy is ex- 
pressed by the probability of pair configurations. 30 
In the present situation the probability of a configura- 
tion on a pair of triangles is approximated as a product 



of t{r\Q, r\\,r\2) variables. The calculation leads to the fol- 
lowing formula: 

^CVAf ^ _g^^ ^ t(^r]o,T]i,r]2)logt{r]o,r]i,r]2) 
+kB ^ Xi,,^ ^ p^,^{riy,rir) logpoi(?7 
+kBYy^Ysy{r]y)logSy{T]y) (24) 



where xqi = a;o2 = J/i = y2 = 2 and X12 =yo = l- 

Substituting Eqs. (|2|) and (|2|) into (|l|) one obtains 
the Gibbs potential as a function of t(r?o, r]i i V2) variables. 
As is well known, Q reaches its minimum at the equilib- 
rium values of these variables. In fact we have seven 
independent variables because the eighth configuration 
probability is determined by the condition of normaliza- 
tion. The numerical treatment is similar to those fol- 
lowed in mean-field approximation. The only difference, 
that instead of the three a^, variables now we have seven 
parameters to minimize the Gibbs potential with respect 
to them. 

All the (mean-field) phase diagrams for R/a — 0.312 
have been recalculated by using this technique. The 
comparisons between the results of mean-field approx- 
imation and CVM are demonstrated in Figs. ^ and 
[1^ . These phase diagrams are typical in the present 
model for the regions 0.1402 < SEt/{e'^/a) < 1.0106 
and SEt > 1.0106e^/a respectively. In both cases the 
CVM suggests lower transition temperatures in compar- 
ison with mean-field results. In general we can say that 
the difference is only a few percent for R/a = 0.312; 
however, it becomes larger for smaller R when the short 
range interactions are more repulsive. 

In brief, this version of CVM confirms qualitatively the 
results of mean-field approximation and leaves its main 
conclusions unchanged. 
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VI. SUMMARY AND CONCLUSIONS 

We have developed a lattice-gas model for investigat- 
ing the formation of different alkali intercalated fuUerides 
with FCC structure. In this (Ising type) model the tetra- 
hedral and octahedral interstitial sites of the cage lattice 
are empty or singly occupied by one type of alkali ion. It 
is assumed that the Coulomb interaction between two al- 
kali ions is screened out by distributing their s electrons 
uniformly on their nearest neighbor Cgo molecules. This 
plausible assumption results in a short range interaction 
which is very convenient for the lattice-gas formalism. 
For ordered structures this model reproduces the elec- 
trostatic energies found by previous authors. [^,0 At the 
same time this approach allows us to study the effect of 
electrostatic forces on the formation of different ordered 
and disordered structures at finite temperatures. 

Besides the electrostatic energy the model takes into 
account the electronic energy of CgQ" ions as well as 
the van der Waals interaction between an alkali ion and 
a Cgo molecule. A series of UHF-MNDO calculations 
shows that the total energy of a charged Ceo molecule 
has a term proportional to the square of its charge. This 
nonlinear contribution mediates an interaction between 
those alkali ions which transfer electrons to the same Cgo 
molecule. This contribution is equivalent to the electro- 
static energy of a charged spherical shell with radius R. 
The extra site energy (SEt) characteristic of the van der 
Waals interaction at the smaller (tetrahedral) interstitial 
voids makes a distinction between the different types of 
alkali atoms. Without this term the tetrahedral sites are 
preferred in the system. In agreement with previous ex- 
pectations this energy contribution increases with 
the size of alkali atom and the octahedral sites are pri- 
marily occupied when 6Et exceeds a threshold value. 

The thermodynamic properties of this model are in- 
vestigated by using a three-sublattice mean-field approx- 
imation. The general features are illustrated via a series 
of phase diagrams. These diagrams draw a picture of 



the effect of van der Waals interaction and the electronic 
energy of the Cg,^ ion on the ordering process. These 
results are confirmed qualitatively by CVM. 

In the present approach the dielectric constant s is 
considered as a parameter. A possible value (e — 9) may 
be estimated by fitting the experimental Raman shift in 
the AgCeo fuUerides. Different values (e « 4) are 
suggested by Pederson and Quong [22| from the polariz- 
ability of Ceo molecules. Desregarding the above uncer- 
tainty the maximum of the predicted ordering tempera- 
tures agree qualitatively with those expected on the basis 
of recent experiments. There are more reliable data on 
the phase separation temperature observed in the K^^Ceo 
system. The present model can describe such type of 
phase diagrams (see Fig. ^ ) in a very narrow range of SEt 
in agreement with the fact that this feature has not been 
observed for other fuUerides. [^3|] Furthermore, applying 
pressure in order to contract the lattice the increasing 6Et 
will couse the sharp decrease of the euthectic tempera- 
ture. Further increase of the pressure can even eliminate 
the phase separation and possibly an orthorombic dis- 
tortion of the cubic structure will appear typical of the 
AiCeo systems at low temperature. 

The ground state energy as a function of SEt and the 
series of phase diagrams describe clearly that the A2Ceo 
composition can appear only for sufficiently low SEt val- 
ues corresponding to small atomic sizes as observed for 
A=Li and Na. 

The generalization of the present description for alka- 
line earth (or other) metals is straigthforward. If the 
intercalated atoms have z valence electrons we should 
substitute a new energy unit, (ze)^/(ea) for e^/(ea) and 
find a suitable SEt expressed in this unit. Due to the 
enhanced energy unit the A2Ceo composition can appear 
in a wider range of SEt- At the same time the stronger 
Coulomb forces contract the lattice and modify the ratio 
R/a which can lead to some change in the phase dia- 
grams. 

The present approach is based on a simplified screening 
mechanism. This hypothesis may be very fruitful for in- 
vestigating other intercalation compounds. The analysis 
of alkali intercalated fuUerides with body-centered-cubic 
and -tetragonal structure is in progress. 
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